Full characterization of Gaussian bipartite entangled states by a single homodyne detector 
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We present the full experimental reconstruction of Gaussian entangled states generated by a type-II optical 
parametric oscillator (OPO) below threshold. Our scheme provides the entire covariance matrix using a single 
homodyne detector and allows for the complete characterization of bipartite Gaussian states, including the 
evaluation of purity, entanglement and nonclassical photon correlations, without a priori assumptions on the 
state under investigation. Our results show that single homodyne schemes are convenient and robust setups 
for the full characterization of OPO signals and represent a tool for quantum technology based on continuous 
variable entanglement. 

PACS numbers: 03.67.Mn, 03.65.Wj, 42.65.Yj 



Introduction — In this letter we address the complete ex- 
perimental characterization of bipartite Gaussian entangled 
states. In our experiment continuous- wave (CW) entangled 
light beams are generated by a single type-II optical para- 
metric oscillator (OPO) below threshold, and then their co- 
variance matrix (CM) is fully reconstructed using a novel 
scheme ^ that involves a single homodyne detector ^ . To 
our knowledge this is the first complete characterization of 
OPO signals without a priori assumptions and paves the way 
to a deeper investigation of continuous variable entanglement 
without experimental loopholes. 

Light beams endowed with nonclassical correlations ISJ are 
crucial resources for quantum technology and find applica- 
tions in quantum communication pi], imaging (5] and preci- 
sion measurement [6, 7] . Their full characterization has a fun- 
damental interest in its own and represents a tool for the de- 
sign of quantum information processing protocols in realistic 
conditions. Remarkably, entangled states produced by OPOs 
are Gaussian states Isl ^ and thus may be fully character- 
ized by the first two statistical moments of the field modes. In 
turn, the CM contains the complete information about entan- 
glement nil , i.e. about their performances as a resource 
for quantum technology. 

Bipartite entangled states may be generated by mixing at a 
beam splitter (BS) two squeezed beams obtained by a degen- 
erate OPO below threshold [12]. The beams exiting the BS 
are entangled [fis!] and a partial reconstruction of the corre- 
sponding CM has been obtained The complete recon- 
struction of a CM has been obtained for different entangled 
states, by varying single-mode squeezing Hs!]. In this config- 
uration two OPOs are used and the amount of entanglement 
critically depends on the symmetry between the two squeezed 
beams. From the experimental point of view this requires an 
accurate setting on the two squeezers and a strict control on 
the relative phase. The measured CM presents some unex- 
pected deviations from a proper form so leaving a question 



open on the reliabilit y of double homodyne schemes due to 
technical difficulties Ill4 A more direct way to gener- 
ate quadrature entanglement is to use a single non-degenerate 
OPO fnll which represents a robust and reliable source of 
EPR-type correlation either below (isl [l^ or above thresh- 
old lEo, [m m. Partial reconstructions of the CM in the 
pulsed regime has been achieved for the spectrally degener- 
ate but spatially non-degenerate twin beams at the output of a 
type-I parametric amplifier 1E3I1 . whereas in the CW regime 
cross polarized beams emitted by a self-locked type-II OPO 
have been examined at frequency degeneracy 1II6I1 . Although 
the correlation properties of OPO signals have been widely in- 
vestigated, no proper CM reconstruction has been performed 
so far. In turn, in previous proposals and experiments non- 
physical CM entries El El [3 12311 or deviations from a 
proper CM [15] appeared, thus requiring a priori hypothesis 
on the measured state to understand the experimental results. 

In this letter we report the first complete measurement of 
the CM for the output of a single non-degenerate OPO. The 
two entangled beams are emitted with orthogonal polarization 
and degenerate frequency by a CW type-II OPO below thresh- 
old. In order to reconstruct the ten independent elements of 
the CM, the beams are optically combined into six auxiliary 
modes, whose quadratures are measured using a single ho- 
modyne detector [if]. The first two moments of the relevant 
quadratures are obtained by tomographic reconstruction us- 
ing the whole homodyne data set, the CM is, then, fully re- 
constructed after assessing the Gaussian character of the sig- 
nal and compared with a general model describing a realistic 
OPO. Entanglement is demonstrated usins the partial trans- 
pose method BIQII . the Duan inequality 111 ill and the stricter 
EPR criterion ll24ll . and quantified upon evaluating the loga- 
rithmic negativity and the entanglement of formation (EoF). 
We also reconstruct the joint photon number distribution and 
demonstrate nonclassical photon correlations by evaluating 
the noise reduction factor. In the following, after defining no- 
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tation and a brief summary of the reconstruction method, we 
describe in details the apparatus and the experimental results. 

Reconstruction method — Upon introducing the vector R = 
(xi, X2, ^2) of canonical operators, in terms of the mode 
operators a/e, Xk = (a^ + ak)/V^, yu = - ak)/V^, 
k = 1, 2, the CM tr of a bipartite state g is defined as the block 
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where A, B and C are 2 x 2 real matrices, (O) =Tt{qO) and 
{fi9} = fg-^gf- In the following we will use the notation 
a = ai and b = a2 and also consider the four additional auxil- 
iary modes c= {a^b)/V2,d= {a-b)/V2,e = {ia^b)/V2, 
and / = {ia — b)/ \pl obtained by the action of polarizing 
beam splitters (PBS) and phase- shifters on modes a and b. 
Positivity of the density matrix for physical states is written 
in terms of the uncertainty relation for the minimum symplec- 
tic eigenvalue V- of the CM, i.e. V- > 1/2. For Gaussian 
states, the state purity is given by fi{a) = {4:y^Det[a])~^ 
whereas separability corresponds to positivity of the partially 
transpose (PPT) density matrix. A bipartite Gaussian state 
is separable iff z>_ > 1/2, where z>_ is the minimum sym- 
plectic eigenvalue of AtrA, A = Diag[l, 1,1,-1]. A con- 
venient measure of entanglement is thus given by the loga- 
rithmic negativity Ef^{cr) = max(0, — ln2z>_ ) 11251] and the 
EoF Ejr{a) can be evaluated following Ref. ll26ll . In addi- 
tion to Ej\f{cr) the Duan criterion gives a necessary condition 
for non- separability in terms of the noise properties of c and d 
I lin whereas a stricter condition [24], referred to as EPR crite- 
rion, involves the conditional variances on Xa and ya obtained 
from a measurement of Xh and y^ for the explicit expressions 
in terms of the noise on modes a, 6, c and d. 

In our experiment, the block A of the CM is retrieved by 
measuring the single-mode quadratures of mode a: the vari- 
ances of Xa and ya give the diagonal elements, while the off 
diagonal ones are obtained from the additional quadratures 

Za {Xa + ya) /V^ and ta = (Xa- ya) / aS (712 

CT21 = ^{{zD - {tl)) - {xa){ya) IM- The block 5 is recon- 
structed in the same way from the quadratures of b, whereas 
the elements of the block C are obtained from the quadra- 
tures of the auxiliary modes c, d, e and / as follows ais = 

{xa){xb),(Ji4 = \{{yl)-{y)))-{xa){yh), 

{xl))-{ya){x,),a2, = \{{yl)-{yl))- 
{ya){yb)- Notice that the measurement of the /-quadratures 
is not mandatory, since {x'j) = {xD -\- (ya)'^ — {xD and 



C^23 = ^{{xj) 



(yj) 



i^a) {y'b) ~ iy'e)' Analogous expressions hold for 
{xl)md {yl). 

In the ideal case the OPO output is in a twin-beam state 
S(C) |0), S(C) = exp{Cat6t - C,ab} being the entangling two- 
mode squeezing operator: the corresponding CM has diago- 
nal blocks A, B, C with the two diagonal elements of each 
block equal in absolute value. In realistic OPOs, cavity and 
crystal losses lead to a mixed state, i.e. to an effective ther- 
mal contribution. In addition, spurious nonlinear processes, 
not perfectly suppressed by the phase matching, may com- 
bine to the down conversion, contributing with local squeez- 



ings. Finally, due to small misalignments of the nonlinear 
crystal, a residual component of the field polarized along a 
may project onto the orthogonal polarization (sa y along 6), 
thus leading to a mixing among the modes OVll . Overall, 
the state at the output is expected to be a zero amplitude 
Gaussian entangled state, whose general form may be written 
as Q, = U(^)S(C)LS(a,6)TLSt(a,6)St(C)Ut(^), 

where T = n r2, with = (1 + nk)~^[nk/{l + n^)]^^^ 
denotes a two-mode thermal state with fik average photons per 
mode, LS(a, 6) = ^(6) ^ ^(6), ^(a) = exp{|(aat2 - 
ikCi^)} denotes local squeezing and U(^) = exp{^a^6 — 
Pab"^} a mixing operator, and 13 being complex num- 

bers. For our configuration, besides a thermal contribution due 
to internal and coupling losses, we expect a relevant entan- 
gling contribution with a small residual local squeezing and, 
as mentioned above, a possible mixing among the modes. The 
CM matrix corresponding to Qg has diagonal blocks A, B, and 
C with possible asymmetries among the diagonal elements. 

Experimental setup — The experimental setup, shown in 
Fig. [II relies on a CW internally frequency doubled Nd:YAG 
laser pumping (@532 nm) a non degenerate OPO based on a 
periodically poled a-cut KTP (PPKTP) crystal (Raicol Crys- 
tals Ltd. on custom design) ll28ll . The use of the a-cut PPKTP 
allows implementing a type-II phase matching with cross po- 
larized signal (a) and idler (b) waves, frequency degenerate 
@1064 nm for a crystal temperature of ^ 53°C. The OPO 
cavity is locked to the pump beam by Pound-Drever tech- 
nique 1I29I1 and adjusted to work in triple resonance by finely 
tuning its geometrical properties 1I27I1 . The cavity output cou- 
pling @ 1064 nm is ^ 0.73, corresponding to an experimental 
line- width of 16 MHz @1064 nm. The measured oscillation 
threshold is Pth ~ 50 mW; during the acquisition the system 
has been operated below threshold at 60% of the threshold 
power. 
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FIG. 1: Experimental setup: A type-II OPO containing a periodi- 
cally poled crystal (PPKTP) is pumped by the second harmonic of 
a Nd:YAG laser. At the OPO output, a half- wave plate (A/2out), a 
quarter- wave plate (A/4out) and a PBSout select the mode for homo- 
dyning. The resulting electronic signal is acquired via a PC module. 

In order to select mode a and b or their combinations c and 
d, the OPO beams are sent to a half- wave plate and a PBS. 
Modes e and / are obtained by inserting an additional quarter- 
wave plate [1]. The PBS output goes to a homodyne detec- 
tor, described in details in jTilioll. exploiting the laser output 
@1064 nm as local oscillator (LO). The overall homodyne 
detection efficiency is 7^ = 0.88 ±0.02. The LO reflects 
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on a piezo-mounted mirror (PZT), which allows varying its 
phase 6. In order to avoid the laser low frequency noise, data 
sampling is moved away from the optical carrier frequency by 
mixing the homodyne current with sinusoidal signal of fre- 
quency 1] = 3 MHz iioll . The resulting current is low-pass 
filtered {B = 300 kHz) and sampled by a PCI acquisition 
board (Gage 14100, IM-points per run, 14 bits resolution). 
The total electronic noise power has been measured to be 
16 dBm below the shot-noise level, corresponding to a sig- 
nal to noise ratio of about 40. 

Reconstruction and experimental results — Acquisition is 
triggered by a linear ramp applied to the PZT and adjusted 
to obtain a 27r variation in 200 ms. Upon spanning the LO 
phase 6>, the quadratures x {0) = xcosO -\- ysinO are mea- 
sured. Calibration with respect to the noise of the vacuum 
state is obtained by acquiring a set of data with the output 
from the OPO obscured. All the expectation values needed 
to reconstruct cr are obtained by quantum tomography (y\\\ . 
which allows to compensate nonunit quantum efficiency and 
to reconstruct any expectation value, including those of spe- 
cific quadratures and their variances, by averaging special pat- 
tern functions over the whole data set. As a preliminary check 
of the procedure, we verified that the CM of the vacuum state 
is consistent with ao = within the experimental errors. 




FIG. 2: (Left): experimental homodyne traces and reconstructed 
CM; (Right): reconstructed Wigner functions. Top plots are for 
modes c and bottom ones for d. Arrows on the homodyne plots 
show the positions of the maximum and minimum variances. 6 is 
the relative phase between the signal and the LO. 

We start our analysis by checking the Gaussian character of 
the OPO signals, i.e. upon evaluating the Kurtosis of homo- 
dyne distribution at fixed phase of the LO ||3Q|i . Besides, we 
checked that the mean values of all the involved quadratures 



are negligible, in agreement with the description of OPO out- 
put as a zero amplitude state. Then, we have measured the 
quadratures of the six modes a-f. We found the modes a and 
b excited in a thermal state, thus confirming the absence of 
relevant local squeezing. Their combinations c, d, e and / 
are squeezed thermal states with squeezing appearing on i/c, 
Xd , te and Zf, respectively. In Fig. [2] we show the experi- 
mental homodyne traces for modes c and d as well as the cor- 
responding Wigner functions, obtained by reconstructing the 
single-mode CM. As it is apparent from the plots both modes 
are squeezed with quadratures noise reduction, corrected for 
nonunit efficiency, of about 2.5 dB. An analogue behavior has 
been observed for modes e and /. The CM of Fig. [2l indeed 
reproduce that of an entangled thermal state with small correc- 
tions due to local squeezing and mixing. The relevant param- 
eters to characterize the corresponding density matrix Qg are 
the mean number of thermal photons fii 0.67, n2 — 0.18 
and entangling photons fls = 2sinh^ \(\ :^ 0.87 ll32n . The 
errors on the CM elements for the blocks A and B are of the 
order Sajk — 0.004 and have been obtained by propagating 
the tomographic errors. In this case phase fluctuations are ir- 
relevant, since the two modes are both excited in a thermal 
state. On the other hand, in evaluating the errors on the el- 
ements of the block C the phase-dependent noise properties 
of the involved modes have to be taken into account, and the 
tomographic error has to be compared with the error due to 
the finite accuracy in setting the LO phase 0. The elements 
(713 and (724 are obtained as combinations of squeezed/anti- 
squeezed variances, which are quite insensitive to fluctuations 
of 6. As a consequence the errors on these elements are given 
by the overall tomographic error 6a jk — 0.004. On the other 
hand, the elements (7i4 and (723 depend on the determination 
of xl J and p which are sensible to phase fluctuations. In 
order to take into account this effect we evaluate errors as the 
fluctuations in the tomographically reconstructed quadratures 
induced by a (56> ^ 20 mrad variation in the LO phase, corre- 
sponding to the experimental phase stability of the homodyne 
detection. The resulting errors are about (5(7i4 = (5(723 — 0.03 
for both CM elements. The off-diagonal elements of the three 
matrices A, B and C are thus zero within their statistical 
errors, in agreement with the expectation for an entangled 
thermal state. As mentioned above, the experimental proce- 
dure may be somehow simplified exploiting the relationships 
among modes, and expressing mode e or / in terms of the 
others: only five modes are then needed. Upon rewriting the 
off-diagonal terms of C in terms of the five modes we arrive 
at (714 = 0.02 ±0.03 and (723 = 0.04 ±0.03 when eliminating 
the mode / and (714 = 0.06 ± 0.03 and (723 = 0.06 ± 0.03 
when eliminating the mode e. Both procedures provide results 
in agreement with those obtained by using the complete set of 
homodyne data for the six modes. 

Since the minimum symplectic eigenvalue of tr is z/_ = 
0.68 ± 0.02 > 0.5, the CM corresponds to a physical state. 
State purity is/i(t7) = 0.31±0.01. The minimum symplectic 
eigenvalue for the partial transpose is z>_ = 0.24 ±0.02, which 
corresponds to a logarithmic negativity Ej\f{a) = 0.73 ±0.02, 
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i.e. the state is entangled, with EoF Ej^{a) = 1.46 ± 0.02. 
In turn, it satisfies the Duan inequahty with the results 0.29 ± 
0.01 < 1/2 and the EPR criterion with 0.21 ± 0.01 < 1/4. 

Entangled Gaussian states as Qg may be endowed with non- 
classical photon number correlations, i.e squeezing in the dif- 
ference photon number. This may be checked upon evaluating 
the noise reduction factor 71 = Ndir{D ah) / {N a + N^) where 
Ndix{Dah) denotes the variance of the difference photocurrent 
Dah = Na — Nh, N = a^a being the number operator, and 
Nk = (Nk) the average photon number. A value 7^ < 1 is 
a marker of nonclassical correlations between the two modes. 
We obtained 7^ = 0.50±0.02, in agreement with the theoret- 
ical description ll33ll for the values of thermal and entangling 
photons reported above. Starting from the CM one can recon- 
struct the full joint photon distribution p(n, m) of the modes 
a and b: the result is shown in Fig. [3] where the correlations 
between the two modes are clearly seen. We have also eval- 
uated the single-mode photon distributions (either from data 
or from the single-mode CM) for modes a-d. Results are re- 
ported in Fig.[3l distributions of a and b are thermal, whereas 
the statistics of modes c and d correctly reproduces the even- 
odd oscillations expected for squeezed thermal states. 




mode a • 
mode c 



p{n) 



|] 



mode b 
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FIG. 3: (Left): Joint photon number distribution p(n, m) for the en- 
tangled state of modes a and b at the output of the OPO. (Right): 
single-mode photon distributions p{n) for modes a and c (top right) 
and b and d (bottom right). The single-mode distributions of modes 
a and b are thermal and correspond to the marginals of p(n, m). The 
distributions for modes c and d are those of squeezed thermal states. 



Conclusion — We have presented the complete reconstruc- 
tion of the CM for the output of a CW type II non-degenerate 
OPO, below threshold and frequency degenerate. The CM el- 
ements have been retrieved as combinations of expectations 
and variances of suitable mode quadratures, obtained by com- 
bining the entangled modes by linear optics. The quantities of 
interest have been obtained tomographically, processing the 
whole data set and thus reducing statistical fluctuations. Upon 
exploiting a general model allowing local squeezing and po- 
larization cross-talking inside the crystal, we have very pre- 
cisely described the experimental CM with the theory under- 
lying parametric downconversion, thus providing a full ex- 
planation of experimental findings. The reconstructed state 
is a Gaussian entangled state close to a two-mode squeezed 



thermal state, the corresponding entanglement and nonclassi- 
cal photon number correlations have been demonstrated. We 
conclude that single homodyne schemes are convenient and 
robust setups for the full characterization of OPO signals and, 
in turn, represent a relevant tool for quantum technology based 
on CV entanglement, e.g., the full characterization of CV 
Gaussian channels by input-output signals' characterization. 
Finally, making use of a single OPO and a single homodyne 
detector, our setup represents also a compact and robust tool 
for entanglement generation and characterization. 
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